Flow behavior of colloidal rod-like viruses in the nematic phase 
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The behavior of a colloidal suspension of rod-like fd viruses in the nematic phase, subjected to 
steady state and transient shear flows is studied. The monodisperse nature of these rods combined 
with relatively small textural contribution to the overall stress make this a suitable model system 
to investigate the effects of flow on the non-equilibrium phase diagram. Transient rheological ex- 
periments are used to determine the critical shear rates at which director tumbling, wagging and 
flow-aligning occurs. The present model system enables us to study the effect of rod concentra- 
tion on these transitions. The results are in quantitatively agreement with the Doi-Edwards-Hess 
model. Moreover, we observe that there is a strong connection between the dynamic transitions and 
structure formation, which is not incorporated in theory. 
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I. INTRODUCTION 

When subjected to shear flow liquid crystals can ex- 
hibit a variety of surprising phenomena, which arise be- 
cause of the anisotropic shape of the constituent rods. 
Theoretically the behavior of a suspension of hard rods 
during shear flow can be described by the equation that 
governs the time development of their probability distri- 
bution function, as derived by Hess and by Doi and 
Edwards [jj. In the absence of a flow, the Doi-Edwards- 
Hess (DEH) theory reduces to the Onsager description 
of equilibrium nematic liquid crystals and can be used to 
describe the isotropic-nematic (I-N) phase transition of a 
hard rod suspension |3l| . The rheological properties are 
predicted to be highly non-linear functions of the Peclct 
number (Pe), which is the ratio of shear rate 7 over ro- 
tational diffusion constant D r . This is not surprising as 
the Pe number can be much larger than unity when the 
rod-like molecules have large aspect ratios. 

The nonlinear response of the rheological properties 
indicates that the shear flow distorts the equilibrium dis- 
tribution of macromolecules or rods. The spatiotemporal 
microstructural changes during flow are even more com- 
plex. At low shear rates, the DEH theory predicts that 
the pseudo vector describing the average alignment of the 
rods, i.e. the "director", undergoes a continuous "tum- 
bling" motion in the plane defined by the velocity and the 
velocity gradient vectors. At high shear rates the director 
is predicted to align with the flow[22|,l2§]. At intermedi- 
ate shear rates it is possible to obtain multiple solutions 
to the Doi-Edwards-Hess equation, which are dependent 
on the initial orientation of the director [lfj, [l2( • For one 
stable solution called "wagging" the nematic director os- 
cillates between two angles in the plane defined by the 
flow and the gradient of the flow. Other solutions such 
as kayaking and log-rolling are also possible, in which the 
director oscillates out of the flow-gradient plane at these 
intermediate shear rates [l5| . 

Experiments on polymeric liquid crystals have con- 
firmed several predictions of the Doi-Edwards equa- 



tion. Using a combination of rheological and rheo-optical 
measurements it was shown that nematic solutions of 
poly(benzyl-glutamate) tumble (PBG) at low shear rate 
and become flow aligning at high shear rates 0. The 
existence of a wagging regime and a potential coexis- 
tence of wagging and log-rolling regimes at intermediate 
flow rates have also been revealed in experiments [l5Ll29| . 
However, there remain significant difficulties when com- 
paring experiments on polymeric liquid crystals (PLC) 
to theoretical predictions. One problem is that different 
levels of the microstructure may lead to different contri- 
butions to the stress tensor . In addition to the molec- 
ular contribution to the stress tensor, textural aspects 
contribute. The latter include Frank elasticity contribu- 
tions due the presence of spatial distortions of nematic 
director, and viscous interactions between 'domains'. In 
addition, there is an indirect effect to the stress tensor 
as the defects disturb the orientation distribution func- 
tion. These textural contributions to the total stress 
dominates the behavior at high concentrations and low 
shear rates [39}, making it difficult to accurately extract 
information about the concentration dependence of dif- 
ferent flow transitions. The textural portion of the stress 
typically displays scaling of the transient rheological re- 
sponse with strain rather than with Pe number |30j . The 
strain scaling is a typical feature of materials where the 
time response is determined by an inherent length scale 
which in the case of PLC's is set by the size of the large 
non-Brownian nematic domains |24|. 

The DEH theory describes the flow behavior of a ho- 
mogeneous ensemble of rods but does not consider any 
polydomain effects. Therefore an ideal system for test- 
ing DEH theory should have small textural contribu- 
tions. In this paper we use rod- like fd virus suspensions 
to access the concentration dependence of the transition 
of tumbling to wagging, and wagging to flow aligning. 
We show that the contribution of textural stress is very 
low, although the spatial distribution of directors still 
has to be accounted for. The main motive for using fd 
virus is the thorough understanding of its equilibrium 
behavior, which has been quantitatively described us- 
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ing the Onsager theory extended to take into account 
the the semi-flexible nature of fd as well as its surface 
charge 32] • Moreover, fd has already successfully been 
used for (micro-)rheology experiments in the isotropic 
phase 0, 0| . The aim of the present paper is to make 
a comparison between the dynamic flow behavior of fd 
suspensions and the available microscopic theoretical pre- 
dictions of the DEH theory for a homogeneous system of 
colloidal rods under shear. 

The paper is organized as follows. In section II wc 
discuss the equation of motion of the orientational dis- 
tribution function and the numerical method we use to 
solve it. The experimental details about sample prepa- 
ration and measurements are given in section III. The 
results are discussed in five parts: the stationary viscos- 
ity of fd suspensions, the concentration and shear rate 
dependence of the oscillatory response to a flow reversal, 
the relaxation after cessation of flow at high concentra- 
tion and in situ microscopy under shear. In section IV 
the textural contribution to the stress tensor is investi- 
gated more detail. Finally we present a non-equilibrium 
phase diagram of shear and concentration dependence of 
different flow transitions. 



II. THEORY 

The distribution of an ensemble of rods can 
be described by the probability density function 
P(ui, ..,iijv, r i, .., r n ) of the positions { r i} and ori- 
entations {Hi} of the rods. Ignoring any spatial cor- 
relations, i.e. restricting to a monodomain, we have 
P{ui, ..,u N , ~ri, .., ~r n ) = pP(ui,..,Uff), where p = 
N/V is the particle density. Therefore, the orientational 
probability density function, or orientational distribution 
function (ODF), fully characterizes the system. The time 
evolution of the ODF for a suspension of rods during flow 
is obtained by solving the equation of motion for the 
ODF, given by the TV-particle Smoluchowski equation. 
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where 1Z(...) — u x V&(...) is the rotation operator with 
respect to the orientation u of a rod. D r is the rotational 
diffusion of a rod at infinite dilution. Furthermore, D is 
the thickness of the rods and L is their length, r = 7 T is 
the velocity-gradient tensor with 7 the shear rate. Here 
we choose 



(2) 




which corresponds to a flow v in the x-direction and its 
gradient Vv in the y-direction. 

The concentration where the isotropic phase becomes 
unstable in the absence of shear flow can be calculated 
by solving the Smoluchowski equation at zero shear rate. 
This equation agrees with Onsagers approach to the I-N 
transition. Often the Maier-Saupe potential is used in- 
stead of the exact potential, which in fact corresponds 
to the first term of the Ginzburg-Landau expansion of 
the outer product in the exact potential given between 
the brackets in Eq. ^ H • Under flow conditions, a rich 
dynamics phase behavior is found as a function of shear 
rate and rod concentration. Marrucci and Maffettone 
were the first to solve the equation of motion of the ODF 
numerically, restricting themselves to two dimensions in 
order to reduce the computational effort [28|. They found 
that the director undergoes a tumbling motion with re- 
spect to the flow direction, resulting in a negative normal 
stress N\ . Larson expanded the ODF in three dimensions 
using spherical harmonics and truncated the expansion 
after checking for convergence j^. This treatment pre- 
dicts a transition from tumbling to " wagging" and finally 
to flow aligning state with increasing shear rates. A clo- 
sure relation is frequently used for the interaction term on 
the right side of Eq. ^ This can greatly bias the results, 
see i.e. Feng et al.[ll|. The location of the flow transi- 
tions in the flow-concentration diagram is very sensitive 
to the choice of the closure, and no satisfactory closure 
has been found up till now. 

In this paper we use a finite element method to nu- 
merically solve the equation of motion for the ODF, thus 
avoiding the use of any specific closure relation. As a 
typical diffusion-convection equation, this equation de- 
scribes the diffusive-convective transport dynamics of an 
orientation of a homogeneous ensemble of thin rigid rods. 
A surface of a sphere is constructed on which a tip of the 
rod moves with respect to its center of mass. The equa- 
tion for the probability of finding the tip of a rod in an 
area is determined by the transport fluxes on its bound- 
aries due to (1) the Brownian diffusion (the first term 
in the brace brackets of Eq. QJ, (2) the convection in- 
duced by the interparticle forces (the second term in the 
brace brackets of Eq. ^| and (3) the convection due to 
the imposed shear flow (the third term of Eq. ^) . 

To solve Eq. n numerically, a discretization scheme is 
used, and meshes on the surface of a unit sphere are con- 
structed. For those operators inside the brace brackets 
which represent the transport fluxes we apply the central 
differences approximations. However, the rotation oper- 
ator outside of the brace brackets needs to be discretized 
using the concept of transport fluxes through the bound- 
aries of the mesh. In other words, the integral form of 
the Eq.(l) is invoked and applied to each of the mesh ele- 
ments. To do this the identity, (u x V&) • F = Vfi • (F x u) 
is used in order to transform the angular transport flux 
of a rod to the translational transport flux of one tip of 
that rod. It differs from the conventional method of dis- 
cretizing a differential equation where the operators are 
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written explicitly into the sum of the first- and the second 
order derivatives and then the latter are approximated by 
selected difference schemes. The advantage of the current 
method is that, since neighboring meshes share bound- 
aries, the fluxes leaving one mesh are always absorbed by 
the surrounding meshes and vice versa. Therefore, there 
is no loss and generation in the total amount of the ODF's 
as the computation proceeds (see Fig^l. In practice a 
40 x 80 mesh was used on the surface of a unit sphere 
with 40 equi-spaced grids in the polar angle and 80 equi- 
spaced grids in the azimuthal angle in a spherical coor- 
dinates. The right hand side of EqQ]is discritized on the 
meshes according to the flux-conservative method men- 
tioned above. A fourth order Adams' predictor-corrector 
method 0] was invoked to follow the time evolution of 
the ODF. More details will be published in a forthcoming 
paper. 

The time-dependent ODF is now used to calculate the 
time-dependence of three parameters characterizing the 
flow behavior of a nematic phase: (1) describing the 
angle between the nematic director and flow direction, 
(2) the scalar magnitude of the director of defined by the 
order parameter Pi and (3) the total stress of an ensemble 
of flowing rods. The angle and magnitude of the order 
parameter are obtained from the order parameter tensor 

S = jdu. uuP(u,t). (3) 

The largest eigenvalue of the order parameter tensor A, 
characterizes the degree of alignment of rods with respect 
to the director given by the corresponding eigenvector n. 
The largest eigenvalue of S is 1/3 in the isotropic phase 
and to 1 for a perfectly aligned nematic phase. Scalar 
order parameter P2 is defined as P2 = (3A — l)/2. 

The stress a 12 is obtained from the deviatoric part of 
the stress tensor derived by Dhont and Briels 

^D=Vo7 + 3pk B T^S-h+^^f + (4) 
ip er [S^:E-iiS:E]|, 

where 

E? = — I du £ du'uu x " X U ' u • u'P(u, t)P(u' t)(5) 
1 3ttJ J |ux u'| 

and 

S (4) = <j> du uuuuP(u,t). (6) 

Here, <f> = jD 2 L~p is the volume fraction of rods, and 
Pe r = j/Dr the rotational Peclet number which is de- 
fined as the shear rate scaled with the rotational diffusion 
of a rod at infinite dilution. The first term between the 
brackets, S — ^1, stems from the Brownian contribution 




FIG. 1: Flux conservation method used in discritizing Eq0 
The rod indicates the orientation of the director with respect 
to the shear flow. The probability of finding a tip of one rod 
in the shaded area of the unit sphere is determined by the 
flux of the probabilities through the boundary of that area. 

to the stress. The second term stems from the direct in- 
teraction between rods and describes the elastic contribu- 
tion to the total stress. The proportionality constant 
is the dimensionless rod concentration and is also called 
the nematic strength. The terms proportional to ~ Pe r 
stem from the flow of the suspension and described the 
viscous contribution to the total stress. This term disap- 
pears instantaneously when the shear is switched off. 

In Fig. |2]we plot the evolution of the three parameters 
(angle 9, order parameter P 2 and stress CT12) as a func- 
tion of strain for different shear rates at a dimensionless 
rod concentration of <f>j^ = 4.5. For this calculation we 
used an initial rod orientation in the flow-gradient plane. 
The flow behavior between Peclet numbers of 4.5 and 
5.0 exhibits a sharp transition from tumbling behavior, 
where the director continuously rotates, to wagging be- 
havior where the director hops back and forth between 
two well defined angles. At higher shear rates the direc- 
tor is found to be flow aligning. The order parameter at 
low shear rates remains unchanged, but is significantly 
reduced at the point of the tumbling to wagging flow 
transition. 



III. MATERIALS AND METHODS 

The viscosity and stress response was measured using 
an ARES strain controlled rheometer (TA instruments, 
Delaware). A double wall Couette geometry was used 
because of fairly low viscosity of the samples. Polarized 
light microscopy images of fd under shear flow were taken 
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FIG. 2: Three plots showing the behavior of the angle of the 
nematic director 9 (a), the magnitude of the nematic order pa- 
rameter^), and the average stress (c) as a function of strain 
after a start up of the flow. The dimensionless rod concen- 
tration is <j)jj = 4.5. Data is obtained by numerically solving 
Eq. fusing the finite element method. The rods are initially 
placed in the flow-gradient plane. For the stress calculation 
only the elastic contribution (Eq. was considered. 



using a Linkam CSS450 plate-plate shear cell. 

The physical characteristics of the bacteriophage fd are 
its length L — 880 nm, diameter D — 6.6 ram, persis- 
tence length of 2200 nm and a charge per unit length of 
around 10 /nm at pH 8.2 01 . When in solution, fd 
exhibits isotropic, cholesteric, and smectic phases with 
increasing concentration [(J [g]. Fd forms a cholesteric 
phase while the DEH theory is valid for nematic struc- 
tures. In practice nematic and cholesteric phase are lo- 
cally almost identical and the free energy difference be- 
tween these phases is very small 0. In this paper we do 
not distinguish between these two phases. The fd virus 



was prepared according to standard biological protocols 
using XLl-Blue strain of E. coli as the host bacteria [33|. 
The standard yields are approximately 50 mg of fd per 
liter of infected bacteria, and virus is typically grown in 
6 liter batches. The virus is purified by repetitive cen- 
trifugation (108 000 g for 5 hours) and re-dispersed in a 
20 mM Tris-HCl buffer at pH 8.2. 



A. fd as a model hard rod system 



The Onsager theory for hard rod dispersions pre- 
dicts a first order phase transition between a disordered, 
isotropic phase and an orientationally ordered, nematic 
phase. Due to hard core athermal interactions consid- 
ered in the Onsager model, the phase diagram is tem- 
perature independent and the rod concentration is the 
only parameter that determines the location of the I-N 
phase transition. The two points spanning the region 
of isotropic-nematic coexistence are called the binodal 
points. The spinodal point is located at a rod concen- 
tration higher then the lower binodal point and is deter- 
mined by the following condition 0-^ = 4. Fd viruses 
are not true hard rods, due to surface charge and limited 
flexibility. As a consequence, their equilibrium phase be- 
havior differs from the ideal hard rod case described by 
Onsager based theory, e.g. DEH. The finite flexibility of 
fd viruses drives the concentration of the binodal points 
to a 30% higher value when compared to equivalent but 
perfectly stiff hard rods. In addition, flexibility also re- 
duces the value of the order parameter of the coexisting 
nematic phase. For fd the order parameter of the co- 
existing nematic is about 0.65 while Onsager theory for 
hard rods in equilibrium predicts the order parameter of 
0.8 [33. The effect of surface charge is to increase the 
effective diameter of the rod D e g and therefore the ex- 
cluded volume interaction between charged rods. As a 
consequence the charge reduces the real concentration of 
the phase transition |35j. 

For the fd suspension used, the binodal point at high 
rod concentration cin occurs at II mg/ml. After tak- 
ing the effects of flexibility and charge into account it 
was shown that the order parameter of the nematic solu- 
tion of fd is quantitatively described by the extensions of 
the Onsager theory to the semi-flexible case |22j. Hence, 
even though fd is flexible and charged, it can be used to 
quantitatively test predictions of the DEH theory. It is, 
however, a very difficult and until now unfulfilled task to 
incorporate charge and flexibility into a non-equilibrium 
equation of motion such as Eq. ^ Therefore in this paper 
we use data from reference |32j to convert the measured 
concentration of fd to the nematic order parameter of the 
sample. After that we compare experiments and theory 
at the same values of the order parameter. 
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FIG. 3: Stationary viscosity as a function of shear rate for 
four different concentration of fd virus at 11.5, 13, 16 and 25 
mg/ml. 
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FIG. 4: The theoretical time averaged viscosity at a dimen- 
sionless concentration of -^0 = 4.5 with (□) and without (■) 
the viscous contribution of the rods, as calculated by solving 
the equation of motion of the ODF for 20 independent initial 
orientations of the director. The lines indicate the transition 
from tumbling to wagging, and from wagging to flow aligning 
as found from Fig. The results are scaled to the experi- 
ment using a typical concentration of 16; mg/ml for p in Eq. 
|U and ZD? f = 20 s" 1 for the rotational diffusion at infinite 
dilution El El- 



IV. RESULTS 

A. Stationary viscosity 

The measurements of a stationary viscosity as a func- 
tion of the shear rate for different fd concentrations are 
shown in Fig. [21 For the lowest concentrations of fd the 
viscosity decreases continuously with shear rate except 
for a small hesitation at a shear rate of 10 s _1 . This hes- 
itation is similar to what is observed for solutions of PBG 
at low concentration in solvent m-cresol [HIIJ. For fd at 
intermediate concentrations, shear thinning becomes less 
pronounced, the hesitation shifts to higher shear rates 
and turns into a local maximum. For the highest fd con- 



centration almost no shear thinning is observed, only a 
pronounced peak in the viscosity. This shear thickening 
behavior has not been previously reported. 

A hesitation in the shear rate dependence of the viscos- 
ity was predicted theoretically by Larson 22]. It was ar- 
gued that the transition from the tumbling regime to the 
wagging regime implies a broadening of the ODF which 
leads to higher dissipative stresses. The broadening of the 
ODF is illustrated in Fig. |2b As can be seen in Fig. ^ 
it is not straightforward that ODF broadening really has 
an effect on the stress. We calculated the time-dependent 
viscosity by numerically solving the equation of motion of 
the ODF for 20 different initial orientations of the direc- 
tor. From the time-dependent ODFs we calculated the 
viscosity using either only the elastic term or both elastic 
and viscous terms. The viscosity is averaged over all 20 
traces and a tumbling period after the transient start up 
flows have died out. The results are scaled to the exper- 
iment using a typical concentration of 16;mg/ml for p 
in eq. 0] and the value of Df n{ = 20 s~\ taken for the 
rotational diffusion at infinite dilution [2(1 [3^] . Fig. 0] 
the stationary viscosity decreases continuously with in- 
creasing shear rate and only shows a hesitation when the 
viscous contribution to the stress is not included. The 
shear rate where this hesitation occurs corresponds with 
the shear rate where the system nematic ordering is sig- 
nificantly reduced and the transition from tumbling to 
wagging takes place, as can be concluded from Fig. |21 
Comparing the model predictions to the experiments it 
should be noted that the experimentally observed fea- 
tures are much more pronounced. Moreover, there is no 
real reason to leave out the viscous contribution although 
it does obscure the behavior we see in the experiment. 
Still, the maximum in the viscosity is interpreted as a 
signature of the transition from tumbling to a wagging 
state. 

There are three observations to keep in mind when con- 
sidering fd in the nematic phase under shear flow, which 
all point to very low stresses in such systems when com- 
pared to polymeric liquid crystals. First, the viscosity of 
fd in the nematic phase is two to three orders of magni- 
tude lower than the viscosity of typical polymeric liquid 
crystals such as poly (benzyl glutamate) (PBG) [3j| , al- 
though the difference in solvent viscosity is only one order 
of magnitude. Second, the range over which the viscos- 
ity of fd suspension varies is more limited with changing 
shear rate and rod concentration: the viscosity lies be- 
tween 70 times the solvent viscosity for low shear rate and 
low rod concentration and 20 times the solvent viscosity 
for high shear rate and rod concentration. Moreover, the 
viscosity as calculated from the equation of motion of 
the ODF is of the same order as the measured viscosity. 
Third, polymer nematics exhibit negative first normal 
stress differences for certain shear rates as was first ob- 
served for PBG solution [Tj|. This is a direct consequence 
of the tumbling of the nematic director. Attempts have 
been made to measure the first normal stress difference 
for nematic fd solutions but due to very low force the 



6 




higher shear rates, the damping again increases (Fig. [SJ;). 
In order to quantitatively characterize the response to a 
nematic to a flow reversal, the data is fitted to a damped 
sinusoidal superimposed onto a asymptotically decaying 
function of the following form: 



v{t) = Vstat ■ { 1 + Ae T d sin(27r — 



if). 



P 



) (l-6-^.)(7) 



This is an empirical choice, but each variable in the fit 
contains important information about the behavior of 
rods in shear flow. Fig. shows the behavior of fit pa- 
rameters as a function of the shear rate at few selected 
concentrations of fd virus. In this figure we indicate with 
vertical dashed lines the shear rates at which the steady 
state viscosity exhibits a local maximum for four differ- 
ent concentrations. Interestingly, these are exactly the 
same shear rates at which the damping constant as 
well as the tumbling period P show a sharp increase. The 
asymptotic constant b, on the contrary, shows a decrease. 
These features disappear for the highest fd concentration. 
Presumably the three regions showing different flow re- 
versal behavior correspond to tumbling, wagging and flow 
aligning regime. This will be discussed in more detail in 
section IV Bl In the next section we first discuss the con- 
centration and shear rate dependence of the tumbling 
period in the regime where rods exhibit tumbling flow 
behavior. 



FIG. 5: The viscosity of the nematic fd solution in a response 
to a flow reversals. The sample is sheared at shear rate +7 
until the viscosity is equilibrated; at time t=0 the shear rate 
is changed to —7. The concentration of fd is kept constant at 
11.5 mg/ml. The data can be nicely fitted to Eq. The fits 
are not shown for clarity. 

signals were too small to be measured. 

B. Flow reversal experiments 

In flow reversal experiments, the sample is first sheared 
at a constant shear rate in one direction until the steady 
state condition is reached. Subsequently, the direction 
of flow is suddenly reversed while keeping the magnitude 
of shear rate constant. Such experiments have been very 
useful in characterizing and understanding the dynamics 
of sheared liquid crystalline polymers [3(j ■ In the present 
work, flow reversal experiments were performed covering 
a wide range of shear rates and fd concentrations. Typ- 
ical flow reversal experiments are depicted in Fig. [5] for 
a fd concentration of 11.5 mg/ml which corresponds to 
c/cin = 1-05. At the lowest shear rates a damped os- 
cillatory response is obtained which decays within few 
oscillations (Fig. [SJi). Increasing the shear rate results 
in a more pronounced oscillatory response, which damps 
out relatively slowly. The oscillatory response in Fig. [SJd 
is most pronounced at a shear rates of 12 s _1 . At even 



C. Tumbling period as a function of shear rate and 
rod concentration 

DEH theory predicts that as the De (or Pe) num- 
ber is increased, that the 'molecular' period of oscilla- 
tion decreases with increasing shear rate in the tumbling 
regime 22]. This feature was never fully explored, since 
in most polymeric liquid crystals it was found that the 
tumbling period was strain scaling, implying that the re- 
sponse overlaps when the period is scaled with the ap- 
plied shear rate and the stress is normalized by its steady 
state value. The strain scaling arises as a consequence of 
the presence of a large, non-Brownian, length scale in the 
sample that determines the time response, even at rela- 
tively high De (or Pe) numbers. This most probably is 
the domain size characterizing the nematic texture. The 
log-log plot of the tumbling period (T=P/7) as a func- 
tion of the shear rate is shown in Fig. Here the data 
are only shown for a low shear rate region which is associ- 
ated with the tumbling region. Strain scaling, if present, 
would give a slope of -1. However, as can be seen in 
the inset of Fig|7| the reciprocal indicating strain scaling 
is only approached and not reached at the highest rod 
concentration studied here. 

The shear rate dependence of the tumbling period is 
compared to the theoretical prediction for the same rod 
concentration as well as the same order parameter, see 
Fig. |HJ The reason for using the order parameter to as- 
sess the theoretical predictions was discussed at length 
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FIG. 6: a) Steady state viscosity as function of the shear 
rate for fd virus at four different concentrations. All viscosity 
curves exhibit shear thinning at low shear rates followed by a 
local maximum in viscosity, b-d) Behavior of the parameters 
obtained from fitting the response of the shear flow reversal 
experiments to Eq. [7| The vertical lines indicate the local 
maximum in viscosity curves. The local maximum in the 
steady state viscosity curve corresponds to maximum of the 
tumbling period P and damping constant Td and minimum of 
asymptotic constat b in the flow reversal. 



in section IIII Al For purposes of comparison the order 
parameter was obtained form x-ray experiments and the 
value of Z)P nf = 20 s _1 . We emphasize that DEH the- 
ory is microscopic and that there are no adjustable pa- 
rameters in the comparison between theory and experi- 
ments. Clearly there is a qualitative correspondence be- 
tween theory and experiment, both showing a continuous 
decrease of the period. The quantitative correspondence, 
on the other hand, is limited. This is probably due to 
fact that texture, although not dominating the response, 
is still present. It will be shown later in section IV Bl that 
the shear rate and rod concentration dependence of a 
tumbling to wagging and wagging to flow-aligning tran- 
sition agree much better with DEH theory. 

The concentration dependence of the tumbling period 
is shown in Fig. [5] Here, theory and experiments are 
compared at a fixed shear rate at which the tumbling 
to wagging flow transition occurs. The tumbling period 
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FIG. 7: The dependence of the tumbling periods on the 
shear rate for different concentrations of the nematic fd. The 
figure shows that the tumbling period scales with a power low 
as a function of the shear rate. The inset shows the power 
law dependence of the tumbling period on the shear rate for 
different fd concentrations. 
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FIG. 8: The period of the oscillations (in units of strain) as a 
function of the Peclet number, where the shear rate is scaled 
with the rotational diffusion of fd at infinite dilution. 



increases with increasing rod concentration (Fig. Et) or, 
equivalently, increasing order parameter of the nematic 
phase (Fig.|5J:)). The increase of the tumbling period with 
increasing order parameter was alread y p redicted using 
a linearized version of the DEH theory[2l| . 

In conclusion, the absence of strain scaling of the tum- 
bling period and the qualitative agreement between the- 
ory and experiment the tumbling period indicates that 
the response of the suspension of /ri-virus is dominated 
by the molecular elasticity arising from the distortion of 
the ODF of particles. 
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FIG. 9: The period of the oscillations (in units of strain) as 
a function of the dimensionless concentration (a) and order 
parameter (b). The shear rate was chosen at the point we 
identify with the tumbling to wagging transition for experi- 
ment (■) and exact theory (o). 
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FIG. 10: The stress relaxation after cessation of flow for 
fd at 25 mg/ml (c/c* = 2.3), varying the initial shear rate. 
The time is scaled by the initial shear rate. The stress is 
normalized by the stress before the cessation of flow. 



cessation of flow, and the time axis is scaled by the shear 
rate. The fast component of the decay takes place at less 
than a tenth of a second, which is comparable to the re- 
sponse of the force re-balanced transducer and therefore 
not shown. The slow component of the stress relaxation 
scales when time is multiplied with the previous shear 
rate, but only from the point that the stress has decayed 
to less than 30 % of its original value, or less for higher 
initial shear rates. From Fig. ^]it can be concluded that 
the contribution to the stress for the highest concentra- 
tion used and for low shear rates is 30 %. This is the 
absolute upper limit for the samples used in this paper. 
It should be noted that for PBG solutions 30 % it was 
found to be the lower limit 391. 



D. Relaxation at high concentration 

In order to measure the relative magnitude of the 
elastic texture contribution to the overall stress, relax- 
ation experiments were performed. For polymeric liquid 
crystals like polybenzylglutamate (PBG) solutions in m- 
cresol, Walker et al. [33] showed that there are three 
different regimes of relaxation behavior, each of which 
is related to a distinct structural relaxation. There is a 
"fast" relaxation of the nematic fluid; a "slower" relax- 
ation that exhibits scaling with the shear rate before the 
cessation of flow, which is due to the indirect contribu- 
tion of the texture to the overall stress; and a "long-time" 
relaxation due to the reorganization of the texture on a 
supra-molecular level which will not be addressed here. 

Stress relaxation experiments were performed in the 
low shear rate "tumbling" region, at shear rates smaller 
than those corresponding to the maximum in the viscos- 
ity. The sample used had a relatively high fd concentra- 
tion of 25 mg/ml, corresponding with c/c/jv = 2.3. Some 
typical responses to the cessation of flow are depicted in 
Fig. The stress is normalized to its value before the 



E. In situ microscopy 

The flow-induced changes of the liquid crystalline tex- 
ture during steady state shear flow were studied using a 
plate-plate geometry in combination with a polarization 
microscope. Measurements were performed for fd con- 
centrations of 14 mg/ml and 25 mg/ml. Typical images 
are shown in Fig. ^2 for different shear rates. Interest- 
ingly the characteristic size of the "domains" was very 
large. Bircfringent regions of up to half a millimeter were 
observed under static conditions. When the sample is 
subjected to shear flow, these domains will elongate and 
eventually disappear, at values of the shear rate which 
correspond to the maximum in the viscosity (see FiglfJ-). 
An important difference between the two concentrations 
is that the elongated domains merge into bands for high 
rod concentration, whereas for the low concentration the 
structure disappears before such bands are formed. In- 
terestingly, this transition to a banded structure in the 
high concentration fluid takes place at a shear rate which 
is higher than the shear rate where the low concentration 
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fluid loses its features. 



V. DISCUSSION 

When comparing the flow behavior of the polymeric 
nematic phase and the colloidal nematic phase of the 
dispersed fd viruses, the most striking observation is the 
qualitative agreement between the two systems, despite 
the fact that fd is an order of magnitude larger. The vis- 
cosity of the fd nematic is much smaller and the rotational 
diffusion of fd is much slower when compared to poly- 
meric liquid crystals. Flow reversal experiments reveal 
typical transitions in the transient rheological behavior: 
damped oscillations occur at low shear rates changing to 
undamped oscillations at intermediate shear rate, which 
disappear if the shear rate is increased even further; the 
time scale of the oscillations of the stress transients is 
comparable. Also other well known phenomena like the 
formation of very large bands upon cessation of flow along 
the vorticity direction which have been studied in detail 
in polymeric systems |37^ , can also be observed here (data 
not shown). 

Having established that fd virus dispersions indeed un- 
dergo a tumbling motion under flow, the dynamic be- 
havior of fd suspensions can be rationalized on the basis 
of the microscopic theoretical predictions for a homoge- 
neous system of rods under shear. Doing so, one impor- 
tant prerequisite needs to be fulfilled, namely that the 
dominating contribution to the stress is coming from the 
nematic fluid and not from the texture. It will be argued 
here that this indeed is the case. Having done so, we 
will be able to map out a phase diagram of the dynamic 
transitions from tumbling to wagging to flow aligning. 



A. Textural evolution during flow 

The word "texture" refers to disclination points and 
lines where the director of the nematic phase changes 
discontinuously, marking domains in the sample. When 
a system containing these domains and disclinations is 
subjected to shear flow, part of the dissipated energy 
is used to destroy these structures. Fig. ^2 shows that 
the domains tend to elongate and align with the flow. 
Disclinations can also cause a direct contribution to the 
total stress resulting in a high viscosity and a very pro- 
nounced shear thinning behavior, typically referred to as 
region I |40|. Experiments on polymeric liquid crystals 
have revealed several features of the flow behavior of ne- 
matic liquid crystals which are attributed to the presence 
of texture in the nematic phase. Tumbling induces dis- 
tortions in the director field and the defects arrest the 
tumbling, thereby inducing an elastic stress. The length 
scale over which this distortion occurs, i.e. the 'domain' 
length scale, is an inherent non-Brownian length scale, 
see ref. 0. As a consequence, stress patterns during flow 
reversal will display strain scaling. Also the damping of 



at rest 
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200 



FIG. 11: The polarization images of the nematic fd at 14 and 
25 mg/ml for a range of different steady state shear rates. The 
dashed line in the bottom right image indicates the border 
between the structured and unstructured regions. 
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the oscillations is explained on the base of the presence 
of the polydomain structure, where e.g. the "friction" 
between the domains would lead to a damping of the 
oscillations The scaling of the stress relaxation 

process after the flow is stopped with shear rate has been 
explained using the same arguments. From such an ex- 
periment the relative contribution to the total stress of a 
homogeneous nematic phase and the polydomain texture 
can be estimated since the relaxation dynamics of the 
nematic phase is much faster than that of polydomain 
structure |39|. 

The micrographs in Fig. II II clearly reveal that texture 
under flow exists in nematic fd dispersions. Their con- 
tribution to the rheology is, however, far less prominent 
when compared to polymeric liquid crystals such as PBG. 
This we can infer from several observations. First, very 
moderate shear thinning is observed in the low shear rate 
regime for the low concentrations, which gradually dis- 
appears with increasing concentration (Fig. Eb)- This 
is very similar to theoretical predictions for a homoge- 
neous nematic phase (Fig. 2b in Ref. |j|). Also, the 
calculated and measured viscosity are of the same order 
of magnitude. In contrast, shear thinning can be fairly 
strong in the low shear rate region (Region I) where tex- 
ture dominates the response and it will increase with 
increasing concentration |27|], although also other mi- 
crostructural features can contribute here [2j. Second, 
the tumbling period is not strain scaling (Fig. [7J, which 
could be due to either a smaller relative magnitude of 
the textural stress or due to the fact that we are not in 
a low enough Pe regime. Third, 'strain' scaling is re- 
covered for the slow 'textural' relaxation process after 
the flow has been stopped. This experiment shows that 
at the highest rod concentrations used and at low shear 
rates the distortional textural contribution is about 30 %. 
For most experiments done this value is probably signifi- 
cantly lower. So, where texture is important, even domi- 
nating the stress response for molecular LCPs, molecular 
elasticity is far more dominating the /d-virus. Though 
we just argued that the texture does not dominate the 
shear response of the system, this does not mean that the 
shear response is not influenced by texture. For one the 
oscillations we observe are still strongly damped, and the 
damping only decreases when the transition to the flow 
aligning state is reached (see the behavior of Td in Fig. HJJ). 
Moreover, the presence of texture might explain the dis- 
crepancy in the behavior of the period of the oscillations 
between experiment and theory (Fig. Most impor- 
tantly, we know from microscopy that texture is present 
under shear (see Fig. Illfl . It should be noted, however, 
that the size of the polydomain structure of the fd dis- 
persions is one order of magnitude bigger as compared 
to PBG [33, so that the density of disclination lines and 
points is about three orders of magnitude lower for fd. 
Note that the length scale of the texture during flow is 
still small compared to the dimension of the flow cell. 
Since the contribution of texture scales with the density 
of the disclinations [27j , texture will be far more dominat- 



ing for e.g. PBG than for fd, even when elastic constants 
are almost the same for the two systems Q,i36]). 



B. Phase diagram of dynamical flow transitions 

In this section the experimental results or combined 
and a non-equilibrium phase diagram of fd rods under 
shear flow is presented. The results for the four fit pa- 
rameters plotted in Fig. show clear transitions at well 
defined shear rates for all fd concentration. Although 
they only give an indirect proof of the transitions, they 
can be used to infer information about the flow tran- 
sitions. For all fd concentrations (except for the high- 
est one) the shear rate where the maximum viscosity is 
reached is identical with the shear rate where the period 
as well as the damping constant start to increase (indi- 
cated by the vertical dashed lines in Fig. [fjj) . The micro- 
scopic observations are in fairly good agreement with the 
transitions inferred from the rheology. Upon approach- 
ing the tumbling to wagging transition from tumbling 
to flow aligning, the texture becomes to faint to resolve 
in the microscope and texture subsequently disappears 
upon reaching the FA region. For the high fd concentra- 
tion, i.e. the sample showing shear banding fFig. llll last). 
one can identify a sharp transition from a structured to 
an unstructured region in the same micrograph. Since 
this picture was taken in the plate-plate geometry, there 
is a shear rate distribution across the image: the shear 
rate is increasing going from the left side to the right. A 
sharp spatial transition therefore also represents a sharp 
transition at a given shear rate. Although, due to the 
method of zero gap-setting, the value of the shear rate 
is not exactly known (±20%), one can still identify the 
shear rate where structure disappears as the shear rate 
where the viscosity reaches its local maximum (the down 
pointing triangles in Fig. HJ-). For low fd concentration 
of (14 mg/ml) the structure disappears around the point 
where the viscosity reaches its local maximum, although 
the morphological transition for the lower concentration 
is less abrupt. 

Fig. E| shows the behavior of flow transitions as a 
function of shear rate for various fd concentrations. For 
the experiment we plotted the Peclet numbers where the 
viscosity shows a local maximum and where the damping 
constant reaches a maximum. The theoretical predictions 
for the tumbling to wagging and wagging to flow aligning 
transitions are obtained from the plots the angle of the 
nematic director 9 under flow, see Fig. [3 Similar to the 
method used in Fig. the experimental concentration 
is scaled to the theoretical concentration in two different 
ways: effective concentration (Fig. ^|a) and the order 
parameter P2 (Fig. ^|b). This figure was shown in a pre- 
liminary paper without a detailed explanation |26j . The 
shear rate is rescaled to the Peclet number by using the 
rotational diffusion coefficient at infinite dilution. Fig. IT2l 
allows us to draw some important conclusions. First, it is 
clear that scaling the concentration with the equilibrium 
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FIG. 12: Phase diagram of flow transitions for the nematic 
fd phase as a function of dimensionless concentration (a) and 
order parameter (b). The experimental points indicate the 
Peclet numbers where the viscosity shows a local maximum 
(■) and where the damping constant Td (•) reaches a max- 
imum. The theoretical points indicate tumbling to wagging 
(□) and wagging to flow aligning (o) transitions. 



order parameter gives better agreement when compared 
to the scaling by the dimensionless concentration. The 
fact that theory and experiment agree without using any 
fitting parameters (P2 was obtained in a separate exper- 
iments [32|) leads to the conclusion that the DEH theory 
describes the flow behavior of the fd nematics quite well, 
as long as the effects of flexibility and charge of the ex- 
perimental rods arc included in the calculation of the 
order parameter. Less convincing agreement is obtained 
when comparing the experimental and theoretically cal- 
culated periods (Fig. |§J). The reason for this could be 
the remaining textural contribution to the overall stress 
which, although small, cannot be neglected. Since we 
deduce from Fig. 112b . that a dimensionless concentra- 



tion of 4>jy — 4 corresponds with a fd concentration of 
16 mg/ml, we used this number a posteriori to scale 
the calculated molecular viscosity in Fig. The pure 
elastic contribution shows a very nice quantitative cor- 
respondence with the experimental data. Interestingly, 
when the viscous term is added, the theoretical viscosity 
is higher than the experimental viscosity, despite of the 
fact that no hydrodynamics is incorporated. 

In the previous subsection it was argued that the in- 
fluence of textural contribution to the stress tensor of fd 
are relatively small, as compared to PLC's. There are 
however strong indications that the dynamic behavior is 
influenced by the macroscopic bands which are formed 
for the samples at the highest concentrations used (see 
Fig. ^]cnd). As can be seen in Fig. the typical fea- 
tures for the transition to wagging disappear: there is no 
increase in the damping constant, nor in the period of the 
oscillations. Moreover, the theory shows only a moderate 
hesitation of the stationary viscosity (which even disap- 
pears then the viscous term is added, Fig. whereas in 
experiments a local peak is observed which is more pro- 
nounced with increasing concentration. The microscopy 
pictures results that at high concentrations the systems 
finds another way to handle the distortion of the particle 
distribution at high shear rates by forming shear bands 
where the overall orientational distribution is alternating, 
as was already observed and partially explained for the 
polymeric systems [23|, |2j| Ijiil] • I n the present work, the 
concentration dependence of the phenomenon at hand 
suggests that this merits further experimental as well as 
theoretical work. In this context one should not forget 
that we compar experiments on charged and semi-flexible 
fd with theory for hard and stiff rods. It could well be 
that these factors also play an important role. It will be 
a major challenge especially to take the semi-flexibility 
into account in the equation of motion. 



VI. CONCLUSIONS 

Colloidal suspensions of rod-like fd viruses are an ideal 
model system to study the behavior of the nematic liquid 
crystalline phase under shear flow. Flow reversal exper- 
iments show signatures for tumbling, wagging, and flow 
aligning behavior, very similar to the behavior found in 
polymeric liquid crystals. The rigid rod nature of the 
fd suspension, possibly combined with a smaller relative 
textural contributions to the overall stress tensor make 
/d-virus a suitable model system for the DEH theory. Im- 
portant in this respect is that the overall viscosity is only 
one to two orders of magnitude higher than the solvent 
viscosity. Also it is important to note that stress relax- 
ation experiments combined with the absence of strain 
scaling in flow reversal experiments suggest that there 
is only a limited contribution of textural aspects to the 
overall stress, even for the highest fd concentration used 
in this work. The shear thickening of the viscosity ob- 
served for a range of fd concentrations is as yet, unex- 
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plained. The maximum in the viscosity occurs at the 
critical shear rate where the tumbling to wagging transi- 
tion takes place. Microscopic observations show that at 
this shear rate the morphological features disappear, sug- 
gesting a strong connection between the dynamic transi- 
tions and structure formation. 

The experimental results have been compared to a mi- 
croscopic theory for rod like molecules subjected to shear 
flow. A non-equilibrium phase diagram is constructed, 
describing the transitions from tumbling to wagging and 
from wagging to flow-aligning as a function of rod concen- 
tration and applied shear stress. When scaling the results 
to the concentration where the isotropic-nematic transi- 
tion takes place, the experiment and theory show only a 
qualitative agreement, possibly due to the fact that the 
real rods are are both semi-flexible and charged. How- 
ever, when scaling the results using the order parameter, 
which is determined by the interactions between the rods, 
theory and experiment show an excellent agreement with- 
out using any fit parameters. Thus, it can be concluded 
that the Doi-Edward-Hess theory accurately captures the 
dynamic features of a hard rod system. Fd dispersions 



constitute such a hard rod system as long as flexibility 
and charge are properly taken into account, which can 
be simply achieved by using the order parameter to scale 
the data. More theoretical work is needed, however, to 
explain the clear connection between the observed band 
formation at high concentrations and the dynamic tran- 
sitions, and to incorporate the effect of flexibility of the 
rods. 
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